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NUMBER VARIANCE OF RANDOM ZEROS 


BERNARD SHIFFMAN AND STEVE ZELDITCH 

Abstract. The main results of this article are asymptotic formulas for the variance of 
the number of zeros of a Gaussian random polynomial of degree N in an open set U C C 
as the degree N —> oo, and more generally for the zeros of random holomorphic sections 
of high powers of any positive line bundle over any Riemann surface. The formulas were 
conjectured in special cases by Forrester and Honner. In higher dimensions, we give similar 
formulas for the variance of the volume inside a domain U of the zero hypersurface of a 
random holomorphic section of a high power of a positive line bundle over any compact 
Kahler manifold. These results generalize the variance asymptotics of Sodin and Tsirelson 
for special model ensembles of chaotic analytic functions in one variable to any ample line 
bundle and Riemann surface. We also combine our methods with those of Sodin-Tsirelson 
to generalize their asymptotic normality results for smoothed number statistics. 
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1. Introduction 

This article is concerned with number and volume statistics for Gaussian random holo¬ 
morphic functions (and sections). To introduce our subject, let us start with the simplest 
case of holomorphic polynomials Pn of degree N of one complex variable. By homogenizing, 
we may identify the space Vn of polynomials of degree N with the space If 0 (CP 1 , O(N)) 
of holomorphic sections of the TV-th power of the hyperplane section bundle over CP 1 . This 
space carries a natural S'U(2)-invariant inner product and associated Gaussian measure 7^. 
To each polynomial p^ we associate its zero set Z PN C CP 1 and thus obtain a random 
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point process on CP 1 . Given an open subset U C CP 1 , we define the integer-valued random 
variable 

■Nn(Pn) = #{- e U : p N (z ) = 0} (1) 

on Vn counting the number of zeros of Pn which lie in U. Clearly, Afjf is discontinuous along 
the set of polynomials having a zero on the boundary dU. It is easy to see from the SU(2) 
invariance that the expected value of this random variable is given by 

e«) = iv ^ a e,, 

where E(X) denotes the expectation of a random variable X and where ©/, is the curvature 
form of the Fubini-Study metric; i.e., the expected zero distribution is uniform on CP 1 with 
respect to its SU( 2) invariant area form. The variance 

Var(A^) = E(A^-E(A ^)) 2 

of A/’v measures the fluctuations of A/’y, i.e. the extent to which the number of zeros of 
individual polynomials conforms to or deviates from the expected number. More generally, 
we can study the same problem for Gaussian random holomorphic sections sjv G H°(M, L N ) 
of powers of any positive holomorphic line bundle L —>■ M over any compact Riemann surface 
M. In [ ISZll , we showed that in this case, the expected value of the random variable A fff has 
the asymptotics 

* E(A ® = sX e ‘ + 0 ® • (2) 

Our first result gives an estimate for the variance of the number of zeros on a domain in 
a compact complex curve, extending and sharpening a result of Forrester and Honner IFHj 
(see also Sodin-Tsirelson Ell): 

Theorem 1.1. Let ( L,h ) be a positive Hermitian holomorphic line bundle over a compact 
complex curve M. We give H°(M, L N ) the Hermitian Gaussian measure induced by h and 
the area form u = |0/j. Let U be a domain in M with piecewise C 2 boundary and no cusps. 
Then for random sections sn G H°(M, L n ), we have 

Length (dU) + 0{N ^ +£ ) . 

This theorem proves a strong form of self-averaging for the number of zeros in U. Here, 
a sequence X N of random variables is called self-averaging if the fluctuations of X N are of 
smaller order than its typical values, or in other words if ^ 0 . 

In higher dimensions, the analogous point process is defined by the simultaneous zeros of 
m polynomials in m variables, or more generally, m sections on an m-dimensional complex 
manifold. In this article we consider instead the simpler ‘volume analogue’ of number statis¬ 
tics for one random polynomial or section s in m dimensions. We let Z s denote the zero set 
of the random holomorphic section s. Recall that the volume of Z s in a domain U is given 
by 

Vol 2m-2[Z s n U}= f 1 uW- 1 = ([Z s ], XU 7 “ pr 

Jz 8 nu(m-1)\ \ (m — 1 ) 



Var(#{z G U : Sjv(z) = 0 }) = VN 
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where [Z s ] denotes the current of integration over Z s . Our higher dimensional generalization 
of Theorem 1 1.11 is the following asymptotic formula for the variance of the volume of the zero 
divisor in a domain with nice boundary. 


Theorem 1.2. Let ( L,h) be a positive Hermitian holomorphic line bundle over a compact 
Kohler manifold (M, to), where uj = We give H°(M, L N ) the Hermitian Gaussian 

measure induced by h, uj (see Definition \2.&) . Let U be a domain in M with piecewise C 2 
boundary and no cusps. Then for random sections sjv G H°(M, L N ), we have 


Var(Vol 2m - 2 [Zs N D U\) = N ~ m+3 / 2 Lvol^!^) + 0(N~^) 


where 


7r m ~ 5/2 

On = —-— C (m- 


Here, we say that U has piecewise C k boundary without cusps if for each boundary point 
Zq G dU, there exists a (not necessarily convex) closed polyhedral cone C C R 2m and a 
C k diffeomorphism p : V —>■ p(V) C M 2m , where V is a neighborhood of Zq , such that 
p{V ft U) = p{V) D C. If M is a complex curve, this condition means that dU is a piecewise 
C k curve with distinct tangents at corners and self-intersection points. 

A model case of Theorem II .21 (as well as of the results stated below) is where M = CP m and 
L = Oil) is the hyperplane section bundle with the SU(m + l)-invariant Hermitian metric h. 
Then sections in H°(M, L N ) are homogeneous degree N holomorphic polynomials on C m+1 , 
and volumes are computed with respect to the Fubini-Study metric oips = fO/i = ^dd log \z \ 2 
on CP m . 

In addition to the number variance problem raised by Forrester-Honner iFHj . the main 
motivation for this article came from the variance and asymptotic normality theorems of 
Sodin-Tsirelson EH for certain model random analytic functions. They consider the smooth 
analogue of number statistics, sometimes called ‘linear statistics’, defined by the random 
variables 

N v(s) = (Z a ,<p)= ( 3 ) 

{z£M:s(z)=0} 

where p> G C^(M) is a test function. In our early paper EZH, we showed that 


E (Z SN , ip) — N f uj A ip + 0(1) , 

J M 

and we gave a crude bound (see EZH Lemma 3.3]) 


( 4 ) 


Var(Zjy, ip) _ ( 1 \ 

on the variance, which was sufficient to prove a strong law of large numbers for the distri¬ 
bution of zeros. In certain model ensembles, Sodin-Tsirelson eh improved this result to a 
sharp estimate as an ingredient in their asymptotic normality result for zeros. Our next re¬ 
sult generalizes their variance asymptotics for the zeros of random polynomials f^ of degree 
N (and their counterparts for model chaotic analytic functions in O(D) and 0(C)) to any 
compact Kahler manifold: 
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Theorem 1.3. Let ( L,h ) be a positive Hermitian holomorphic line bundle over a compact 
Kahler manifold ( M,to ), where u> = 10^. Let (p be a real (2m — 2 )-form on M with C 3 
coefficients. We give H°(M, L N ) the Hermitian Gaussian measure induced by h and the area 
form to (see Definition s. &) . Then for random sections sn € H°(M, L N ), we have 


Var (Z SN ,v)=N~ m 


K m \\ddip\\l + Q(N 2 +e ) 


where 


him 


7r 


m —2 


4 


C(m + 2) . 


Here, ||<9<9</3|| 2 denotes the £ 2 norm of ddtp, i.e. writing idd(p = to m , we have ||<9<9<^||| = 
f = f ifjddip. (Of course, we may assume that ip is of bidegree (m — 1, m — 1), since 

( Z SN , ip) = 0 for forms <p of other bidegrees.) 

In particular, for the case dim M = 1, we note that \dd(p\ = ||Ay?|, and thus 


Var(Z Sjv ,^) =iV 1 


03 ) 

167T 


||A^|| l + 0(N~-^) 


( 6 ) 


The leading term in © was obtained by Soclin and Tsirelson 1ST] for the case of random 
polynomials s^ £ i/°(CP 1 , 0(N )) and random holomorphic functions on C and on the disk. 
(The constant ^0 was given in a private communication from M. Sodin.) 

Our final result is an asymptotic normality result of the type proved in EH. It follows 
very easily from the analysis underlying Theorem 11.31 together with a general asymptotic 
normality result of Sodin-Tsirelson. 


Theorem 1.4. With the same notation and hypotheses as in Theorem } 1 .,4 the distributions 
of the random variables 

(Z SN ,y) - E (Z SN ,y) 
y/Var (Z SN ,<p) 

converge weakly to the standard Gaussian distribution A/”(0,1) as N —> oo. 

We let A/”(0, a) denote the (real) Gaussian distribution of mean zero and variance a 2 . 
Substituting the values of the expectation and variance of (Z SN ,<p) from © and Theorem 
11.31 respectively, we have 


Corollary 1.5. With the same hypotheses as in Theorem, II.,4 the distributions of the 
random variables N m ^ 2 (Z SN — Nto,ip) converge weakly to A/"(0, y/nff 1111 2 ) as N —► oo. 


Let us briefly mention some key ideas in the proofs and the relation of the Sodin-Tsirelson 
methods to ours. The Sodin-Tsirelson estimate was based on their construction of a ‘bi- 
potential’ for the pair correlation measures, i.e. functions Gn(z,w ) such that 

A z A w G n (z, w) = I <2 (z, w). (7) 


Here, Kf is the ‘pair correlation function’ for the zeros of sections of Sn, that is, the 
probability density that a section in Sn has zeros at two points z and w of CP 1 . The bi¬ 
potential is given in (STj as a power series in the Szego kernel for 0(N) —> CP 1 . (Here, the 
notation in © is taken from [BfiZlj and is not used in EDO In fact, the same bipotential 
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already arose in [Rj-jZ 1 j in the setting of line bundles over a compact Kahler manifold as the 
bi-potential for the ‘pair correlation current,’ i.e. 

-d z 3 z d w 8 w G N (z,w) = E ( Z SN (z ) ® Z SN (w )) . (8) 

We build on our analysis of this bi-potential in [BSZ1] to prove Theorem 11.81 

The other main ingredient in the proofs are estimates derived from the off-diagonal asymp¬ 
totics of the Szego kernel in jSZ2j . For the sake of completeness, we review the derivation of 
these off-diagonal asymptotics in 10] Off-diagonal estimates of the Szego kernel with sharper 
(exponentially small) remainder estimates are given in IDLMl IMMj . but the estimates of 
|SZ2| already suffice for our applications. 

Although we are emphasizing positive line bundles over compact Kahler manifolds, our 
results (and their proofs) extend with no essential change to positive line bundles over non¬ 
compact Kahler manifolds for which the orthogonal projector Iljv onto £ 2 H°(M, L N ), the 
L 2 holomorphic sections of a positive line bundle with respect to a Hermitian metric and the 
Kahler volume form, has analytic properties similar to those in the compact case. A model 
for a positive line bundle over a noncompact Kahler manifold is provided by the Heisenberg 
line bundle L H —> C m associated to the reduced Heisenberg group by the identity character, 
as described in detail in [ BSZ2 , IBSZ3| . In this case, the analogue of Theorem 11.21 is an 
asymptotic formula (Corollary 12.511 for the volume variance of the zeros of random Gauss¬ 
ian entire functions on the dilates VN U of a domain U C C m . Other model examples 
are given by homogeneous Hermitian line bundles over bounded symmetric domains with 
curvature equal to the Bergman Kahler metric. We briefly discuss the extension to random 
holomorphic sections in the noncompact case in H2.2I 

It will readily be recognized that the variance and normality problems in higher dimensions 
make sense for the simultaneous zeros of k independent sections si,..., Sk of a line bundle over 
an m-dimensional complex manifold and are perhaps most interesting for the full codimension 
case k = m. The same problem may be posed for the critical points of a single Gaussian 
random section. However, new technical ideas seem to be necessary to obtain limit formula 
for the intersections of the random zero currents Z Sj . We hope to return to this problem 
elsewhere. 

In conclusion, we thank M. Sodin for discussions of his work with B. Tsirelson on number 
variance and asymptotic normality for random analytic functions of one variable. 

2. Expected distribution of zeros and Szego kernels 

In this section, we review the basic formula from BSZ1, BSZ21I3ZT] for the expected 
distribution of zeros of Gaussian random sections of holomorphic line bundles. We state it 
here in a general framework which we shall use in our forthcoming paper on zeros of random 
fewnomials SZ3 j. 

We let (L, h) be a Hermitian holomorphic line bundle over a complex manifold M (not 
necessarily compact), and let S be a finite-dimensional subspace of H°(M,L). We suppose 
that dim S > 2 and we give S a Hermitian inner product. The inner product induces the 
complex Gaussian probability measure 



j= 1 


( 9 ) 
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on where {Sj} is an orthonormal basis for S and dc is 2n-dimensional Lebesgue measure. 
This Gaussian is characterized by the property that the 2 n real variables Recj,Imcj (j = 
1,... ,n) are independent complex Gaussian random variables with mean 0 and variance 1; 
i.e., 

E Cj = 0, E CjCk = 0, E CjCk = Sj k ■ 

We let 

n 

n s (z,z) = £||S,(*)||J, ««, (10) 

3 =1 

denote the Szego kernel for S on the diagonal. (See H2.ll for a discussion of the Szego 
kernel.) We now consider a local holomorphic frame ex over a trivializing chart U, and we 
write Sj = fjCL over U. Any section s G S may then be written as 

n 

s = (c,F)ef N , where F = (/i,..., f k ) , (c, F) — ^ Cjfj . 

3 = 1 

If s — feL, its Hermitian norm is given by \\s(z)\\h = a(z)~^\f(z)\ where 

a (z) = \\e L (z)\\~ h 2 . (11) 

Recall that the curvature form of (L, h) is given locally by 

O h = dd log a , 


and the Chern form Ci(L, h ) is given by 


ci{L,h ) 


v^T 

2vr 


0 


h 


v^T 

2vr 


dd log a . 


( 12 ) 


The current of integration Z s over the zeros of s = (c,F)cl is then given locally by the 
Poincare-Lelong formula: 

Z s = ^-^-dd log | (c, F) | . (13) 

7r 

It is of course independent of the choice of local frame ex and basis {Sj}. 

We now state our formula for the expected zero divisor for the linear system S: 


Proposition 2.1. Let (L, h) be a Hermitian line bundle on a complex manifold M, and let 
S be a finite dimensional sub space of H°(M, L). We give S a Hermitian inner product and 
we let 7 be the induced Gaussian probability measure on S. Then the expected zero current 
of a random section s G S is given by 


E 7 (Z S ) = — — ddlogTls(z,z) + d(L,h) . 

ZTT 

Proof. Let {Sj} be an orthonormal basis of S. As above, we choose a local nonvanishing 
section ex of L over U C M, and we write 

n 

s = ^CjSj = (c, F)e L , 

3 = 1 
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where Sj = fjeL, F = {fi, ■ ■ ■, fk)- As in the proof of m, Proposition 3.1, we then write 
F(x) = \F(x)\u(x) so that |w| = 1 and 

log I (c, F) I = log |F| + log I (c, u) I . (14) 

A key point is that E( log |(c, w)|) is independent of z (and in fact, is a universal constant 
depending only on n) , and hence E(dlog |(c, «)l) = o. 

Thus by (USD, we have 

(E 7 (Z s ),cp) = [ (log | (c, F) |, dd<p) d'y(c) 

77 Jc n 

= f (logIFI,<9(9v?) dy(c) + f (log \{c,u)\,dd(p)d'y(c) , 

77 Jc n 71 J C n 

for all test forms ip G D m_1,m_1 (f/). The first term is independent of c so we may remove 
the Gaussian integral. The vanishing of the second term follows by noting that 



(log I (c, u)\, ddip) d'y(c) 


d'y(c) / log\(c,u)\ddip 

J M 


M L/C” 


log|(c,u)|d 7 (c) 


dd<p = 0 , 


since f log |(c, u) \ d^(c) = ^ J c log|ci|e \ ci][2 dc\ is constant, by the U (ri)-invariance of d'y. 
Fubini’s Theorem can be applied above since 



|log | <c, w) | ddip | dy(c) = 


1 / 

77 J c 


log |Ci 11 e ' Cl ' 2 dc] 



< Too • 


Thus 

F^s) = ^^ d d\og\F\ 2 = ^^ d d(\og S ^\\ s j \\l + \°ga) . (15) 

Recalling that Hs(z,z) = Y ll*Sy (-^) II*. an d that ci(L,h) = (99 log a, the formula of the 

proposition follows. □ 


Remark: The complex manifold M, the line bundle L and space S as well as its inner product 
in Proposition 12 .II are all completely arbitrary. We do not assume that M is compact or that 
(L, h ) has positive curvature. We do not even assume that S is base point free. If S has no 
base points (points where all sections in S vanish), then we have the alternate formula (see 

BZ2) 


E -y(Z s ) — ‘hjo/ps , 


where $5 : M —> P S* is the Kodaira map and a/ps is the Fubini-Study form on P5*. In the 
general case where there are base points, we have 


E 7 (Z S ) — ‘FJa/ps T D , 
where D is the fixed component of the linear system P«S. 
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2.1. Powers of an ample line bundle. We now let L —> M be an ample line bundle on a 
compact complex manifold M. We consider tensor powers L A = L® N of the line bundle, and 
we let S = H°(M, L N ). We further choose a Hermitian metric h on L with strictly positive 
curvature and we give M the Kahler form u = |@h = 7 TCi(L, h). 

We now describe the natural Gaussian probability measures on H°(M, L N ) used in jSZll 
1SZ21 iB^Zll IBSZ2] . For the case of polynomials in one variable, these Gaussian ensembles 
are 


equivalent to the SU(2) ensembles studied in EEL Hanl IN VI SZlj and elsewhere. 


Definition 2.2. Let ( L,h ) —> (M,u) be as above, and let h N denote the Hermitian metric 
on L n induced by h. We give H°(M, L N ) the inner product induced by the Kahler form u ; 
and the Hermitian metric h N : 

(si,s 2 ) — f h N (si,s 2 )—:ca m , Sl ,s 2 eH°(M,L N ) . (16) 

J m ml 

The Hermitian Gaussian measure on H°(M, L N ) is the complex Gaussian probability mea¬ 
sure 7 at induced by the inner product m-- 


d N 

rf 7 w(s) = — e" |c|2 dc, S = Y1 c i S ? , 

77 3 = 1 

where {S± ,..., Sjf N } is an orthonormal basis for H°(M, L N ). 

As in jSZll iB^Zll and elsewhere, we analyze the Szegb kernel for H°(M, L N ) by lifting 
it to the circle bundle X^M of unit vectors in the dual bundle L~ l —> M with respect to 
h. In the standard way (loc. cit.), sections of L N lift to equivariant functions on A". Then 
s G H°(M, L N ) lifts to a CR holomorphic functions on X satisfying s(e l6 x) = e tNe s(x). 
We denote the space of such functions by TL 2 N (X). The Szego projector is the orthogonal 
projector : £ 2 (X) —> TL 2 N {X ), which is given by the Szego kernel 


H N (x,„) = Y.H l - x )Sj'(y) (x.yeX). 

3 = 1 

(Here, the functions S^ are the lifts to Tt 2 N (X) of the orthonormal sections S^ 1 ; they provide 
an orthonormal basis for T( 2 n (X ).) 

Further, the covariant derivative Vs of a section s lifts to the horizontal derivative V^s of 
its equivariant lift s to A"; the horizontal derivative is of the form 



For further discussion and details on lifting sections, we refer to 
We shall write 


SzTj. 


(17) 


|n iV (^, td)| := |n A r(x,r/)| , z = tt(x), w = vr (y) G M . 

In particular, on the diagonal we have Hat(^, z) = where tt(x) = z. Note that 

nyv(- 2 , z) = n 5 ( 0 , z) as defined by dTOT) with S = H°(M, L N ). 
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2.2. Random functions and Szego kernels on noncompact domains. As mentioned 
in the introduction, Theorems I1.1H1.4I extend with no essential change to positive line bun¬ 
dles over noncompact complete Kahler manifolds as long as the orthogonal projection onto 
the space C 2 H°(M, L N ) of L 2 holomorphic sections with respect to the inner product dTHl) 
possesses the analytical properties stated in Theorem 14.11 land mostly proved in JSZ2j ) for 
Szego kernels in the compact case, ft would take us too far afield to discuss in detail the 
properties of Szego kernels and random holomorphic sections in the noncompact setting, but 
we can illustrate the ideas with homogeneous models. 

Before discussing our specific noncompact models, we first note that Proposition 12 .II holds 
for infinite-dimensional spaces of Gaussian random holomorphic sections. There are sev¬ 
eral equivalent ways to describe Gaussian random analytic sections or functions in an in¬ 
finite dimensional space (e.g., {Jaj G.j ISTj ). To take a simple approach, we suppose that 
{Si, S 2 ,...,S n ,...}isa,n infinite sequence of holomorphic sections of a Hermitian line bundle 
(L, h ) on a (noncompact) complex manifold M such that 


sup^HSj 


< +oo 


zeK 


3 =1 


for all compact K C M . 


(18) 


We then consider the ensemble («S, d'y) of sections of L of the form 


S 


oo 

s = > CjSj : c,- G C 


j= 1 




(19) 


i.e. we consider random sections s = fffJLiCjSj, where the c 3 are i.i.d. standard complex 
Gaussian random variables. It is well known that m implies that the series in irrai i almost 
surely converges uniformly on compact sets (see e.g. Jal lKahn . and hence with probability 
one, s E H°(M,L). We then have: 


Proposition 2.3. The expected zero current of the random section s e S in (CHI) is given 
by 


E (Z s ) = G—dSlogn (z,z)+ Cl (L,h) , 

Z7T 

where 

oo 

n(z,z) = £||s,v)||y 

3 = 1 

Proof. The proof is exactly the same as the proof of Proposition 12. II using the ensemble (trail 
with the infinite product measure, except we cannot use unitary invariance to show that 

J log | (c,«) d,~j(c) = i / c log|ci|e- |ci|S dci, /1 log |<c, ti) 11d-y(c) = j f c | log |ci|| e" |0l|2 (ici. 

( 20 ) 

To verify (ESI) in this case, we note that (c,u(z)) is a complex Gaussian random variable of 
mean 0 and variance 1 (see pal EaE]), and hence 


f({c,u))d'y(c) = - / /(C) e d( , 


7T 


The identities of m then follow by letting /(C) = 


for all f e £ 1 (C,e~M 2 d() . 
log |CI, resp. /(C) = I log |C| 


□ 
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We are interested in the case where (L, h) has positive curvature, M is complete with 
respect to the Kahler metric a; = |©/i, and {Sj } is an orthonormal basis of jC 2 H°(M,L) 
with respect to the inner product (USD- Note that with probability one, a random section 
s is not an C 2 section (since ||s||2 = || c||2 = +00 a.s.), but is a holomorphic section of L. 
(Equivalently, C 2 H°(M, L) carries a Gaussian measure in the sense of Bochner-Minlos; see 

KUO ' 

The first model noncompact case is known as the Bargmann-Fock space 


T :=H 2 (C m : e-^ 2 ) = \f C m ) : / |/|V |z|2 dz < +00 


We can regard elements of T as C 2 sections of the trivial bundle Lh over C m with metric 
h = e _ + . The associated circle bundle A" can be identified with the reduced Heisenberg 
group; see [HSzHl §2.3] or [B?Z 2 l §1.3.2], Then 3= = C 2 H°(C m ,L n ) = Hj(X), and more 
generally, 


C 2 H°(C m ,L£) = H 2 N ( X) 


if e 0 { C m ) : 



\ 2 e- N W 2 dz < +00 


( 21 ) 


In dimension one, this example is referred to as the ‘flat model’ in (HD. 

An orthonormal basis for the Hilbert space C 2 H° (C m , L^f) with inner product (/i,/ 2 ) = 
J C n , fiffe~^ 2 dz is 



N m/2 (7Vl/ 2 z) fe ] 

^ m/2 Vk\ J fceNm ’ 


where we use the usual conventions z k = z'l 1 ■ ■ ■ zk\ — k\\ ■ ■ ■ k m \, \k\ = k\ + ■ —h k m . A 
Gaussian random section is defined by 


In(z) = c ^ S k (-) 
fceN m 


]ym/2 

yj-m/2 


E 

ke N m 


Ck 

71 


(N 1 / 2 z) 


k 

? 


( 22 ) 


where the coefficients Ck are independent standard complex Gaussian random variables as in 
©. As mentioned above, the random sections /at are almost surely not in C 2 H°(C m , Lff). 
However, they are almost surely entire functions of finite order 2 in the sense of Nevanlinna 
theory. Indeed, we easily see from © that 


7 ({c G C°° : |c fc | 2 < 2 log kj for kj > 2 j j > 0 , 

and hence it follows from the zero-one law that || 2 = O log kj'j a.s. Therefore, by 
Cauchy-S chwartz, 


l/WI = o(E( 1” £ ) i * i ) 1/2 =o( e <H«»p) 

for all e > 0. Thus we have an upper bound for the Nevanlinna growth function, 


a.s. 


T(f N ,r) ■= Ave {N=r} log + \f(z)\ < sup log + \f(z)\ < 

\z\<r 


N 


+ o(l) 


a.s. 
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(where o(l) denotes a term that goes to 0 as r —> oo, for each fixed N > 1). On the other 
hand, if T(f N ,r ) = 0(e^ _£ ) iv l 2 l 2 ), then f N <E £ 2 //°(C m , Ljj), which has probability zero. 
Thus, 


lim sup 


T(f N ,r) 


N 

~2 


a.s. 


To use the proofs in §(JSHZIto show that Theorems II . lllTD hold for the line bundle Ln, we 
need only to verify that the Szego kernel 11^, i.e. the kernel of the orthogonal projection to 
C 2 H°(C m , Ln), satisfies the diagonal and off-diagonal asymptotics in Theorem 14.11 In the 
model Heisenberg case, the Szego kernel is given by 


n *(z,0-,w,ip) = e iN « > -rt S k (z)S k (w ) 

fceN m 


i\rm 

IV c iN(e-y)+Nz-w-%-(\z\ 2 + \w\' 2 ) 

7 T m 


(23) 


(see IBSZ21) and visibly has these properties. 

Another class of homogeneous examples are the bounded symmetric domains fl C C m , 
equipped with their Bergman metrics u = \dd log K(z, z) where K(z, z ) denotes the Bergman 
kernel function of fh Let (L, h ) —> be the holomorphic homogeneous Hermitian line bundle 
over with curvature (1,1) form u. It was observed by Berezin [Berl that the Szego kernels 
Il N for jC 2 H°(Q, L n ) also have the form C^e Nl ^ z ’ w ' i where Cn is a normalizing constant and 
"0 = log K(z,w). In the case of the unit disc D C C with its Bergman (hyperbolic metric) 
~^dd log(l — |^| 2 ), the space C 2 H°(D, L N ) may be identihed with the holomorphic discrete 
series irreducible representation of SU( 1,1) (cf. |Knl p. 40]), that is with the space of 
holomorphic functions on D with inner product 


\\f\\ 2 N= [ \f{z)\ 2 (l-\z\ 2 ) N ~ 2 dz. 

J D 

The factor comes from the Hermitian metric. An orthonormal basis for the 

holomorphic sections of L N is then given by the monomials (] Ar +”” 1 ]j 1//2 ^ n (n — 0,1, 2,...). 
The Szego kernels are given by njv(z, w) = (1 — zw) N . The Szego kernels also visibly have 
the properties stated in Theorem 14.11 These ensembles are called the hyperbolic model in 
EH- Random SU(1,1) polynomials are studied in |BR| . where further details can be found. 
Thus our proofs also yield the following result: 


Theorem 2.4. Theorems \l.ll\l. 4 \ hold for the zeros of sections in the following ensembles: 

• Gaussian random sections f^ G H°( C m , Lji) given by (1^1) : 

• Gaussian random sections of the holomorphic homogeneous Hermitian line bundle 
(L, h) over a bounded symmetric domain Gl C C m , as described above. 


We note that taking the iV-th power of the line bundle Ln —» C m (i.e., taking the N -th 
power of the metric e~^ z ^ ) corresponds to dilating C m by \fN. Precisely, the map 

t n : C 2 H°(C m ,L u ) - C 2 H\C m ,L£) , (r N f)(z) := N m ' 2 f (TV 1 / 2 z ) , 

is unitary. Thus we can restate our result on the volume (or number, in dimension 1) variance 
for the Bargmann-Fock ensemble as follows: 


/M = £ 

fee N m 


Ck „k 

Vk\^ 


) 


Corollary 2.5. Let 
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where the coefficients Ck are independent complex Gaussian random variables with mean 0 
and variance 1. Let U be a domain in C m with piecewise C 2 boundary and no cusps, and 
consider its dilates Un '■= vNU. Then, 

Var(Vol 2 m _ 2 [Z/ fl Un]) =^ m Vol 2m -i(dU N ) + 0(N~* +e ) , 
where u m = LUJfl £( m + I). 

Note that 8Un = N m ~ l / 2 dU , so we have Var(Vol 2 m _ 2 [Z/ fl y/N U}) ~ A m_1 / 2 . 

Off-diagonal estimates for general Bergman kernels of positive line bundles over complete 
Kaliler manifolds are proved in [MM] using heat kernel methods. The relevant issue for 
this article is the approximation of the L 2 Szego kernel by its Boutet de Monvel -Sjostrand 
parametrix in the noncompact case. The analysis of Szego kernels on noncompact spaces lies 
outside the scope of this article, so we do not state the general results here. But it appears 
that the general results of | IMM: give sufficient control over Szego kernels in the noncompact 
case to allow Theorems 00 to be extended to all positive line bundles over complete 
Kahler manifolds. 


3. A BIPOTENTIAL FOR THE VARIANCE 


Our proofs of Theorems 11 . 11 - ITvil are based on a bipotential implicitly given in 
describe our bipotential Qn{z,w), we define the function 


szn . To 


G(t) := f lQg(1 S) ds , 0 < t < 1 . (24) 

4tt 2 J o s 

Alternately, 

~ i r°° 

G(e~ x ) = - 5-2 / log(l -e~ 2s )ds , A> 0 . (25) 

A 

(The function G is a modification of the function G defined in EZB; see (l4Tll .l We also 
introduce the normalized Szego kernel 


Pn(z,w) 


|rijv(z, w) | 

n Ar (^,z)^n Ar (w;,w)2 


(26) 


Definition 3.1. Let ( L,h) — > (M,cu) be as in Theorem, s \l. L 3U., b A The variance bipotential 
is the function Qn : M X M —> [ 0, + 00 ) given by 

Q N (z , w) = G(P n (z , in)) = Z ds . (27) 

471 "= J 0 s 

We remark here that Qn is C°° off the diagonal, but is only C 1 and not C 2 at all points on 
the diagonal in M x M, as the computations in show. 

The variance in Theorems OO can be given by a double integral of the bipotential, as 
stated in Propositions 13.21 and 13.51 below: 


Proposition 3.2. Let (L, h) —> (M, cn) and ip be a (2m — 2) -form on M with C 2 coefficients. 
Then 



Qn(z,w) (iddip(z)) ( iddp>(w )) 


m J M 


Var(Z,„,*j) 


(28) 
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To begin the proof of the proposition, we write 

= (Sf,.... Sit) 6 »•(«, lN )“" . (29) 

where {S^} is an orthonormal basis of H°(M, L N ). As in the proof of Proposition 12.11 we 
write 

^n(z) = \^f N (z)\u N (z) , 

where |\Ew| := ||5'j V ||feiv) 1/2 , so that \u N \ = 1 . 

We hrst establish a less explicit variance formula: 


Lemma 3.3. 

Var (Z SN ,ip) = ^f f (idd(p(z))(idd(p(w)) f log \{u N (z), c)| log \{u N (w), c)| dr/ N (c) . 

11 JM JM JC d N 

Proof. We write sections % 6 H°(M, L N ) as 


d N 


sn = c i S f = ^ ’ C = (Ci, . . . , c dN ) . 

3 = 1 

Writing = Fef N ", where is a local nonvanishing section of L, and recalling that 


(30) 


u = -@h = -idd log ||ez,|U , 


we 


have by m, 


Zs N = -99 log I (c, F) I = -99log I {c,V N ) I - -99log Hef^Hh 

71 71 71 

i tv 

= -99log |(c, Tat)! H- u. 

7T 7 T 

Let 99 be a test form, and consider the random variable 

Y N := (-ddlog\(c,ty N ) \ , (p) = f log | (c, tEw) | , —ddcp] 


71 


71 J 


so that (Z SN , <f ) — Y n + NC, where C — f M uj A y?; hence 

Var(Z Sjv ,< / 9) = Var(yjv) . 

By (USD, we have 


E(Y/v) = 99 log | Tjv | ,</?=-/ log |^jv| dd<p , 


71 


7T 


'M 


(31) 


(32) 


(33) 


whereas by (E2D, we have 

E ( y n) = [ [ (iddip(z)) (idd<p(w)) [ log |(c, tfjv(z))| log |(c, 'L iV (w))| drm(c) . (34) 

71 Jm iif 

Recalling that = I'EwImn with |wjv| = 1, we have 
log|(T iV (^),c)| log|(^iv(w),c)| = log|^jv(z)| log|^iv(w)|+ log|^ 7 v(z)| log|(ujv(«0,c)| 

+ log|tfjv(w)| log | ('itjv(^), c) | 

+ log I(tijv(w), c>I log|(ttjv(z),c)| , 


( 35 ) 
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which decomposes (EH into four terms. By 4HHD, the first term contributes 


(idd<p(z))(iddip(w))\og\^ N (z)\ log |'I' JV (w)| = (E Y N ) 



11 J M J M 


( 36 ) 


The c-integral in the second term is independent of w and hence the second term vanishes. 
The third term likewise vanishes. Therefore, the fourth term gives the variance Var (Z SN , (p). 

□ 

We now complete the proof of Proposition 13.21 bv evaluating the c-integral of Lemma I’d.Ml 
Lemma 3.4. We have: 

4 / log\(u N (z),c)\log\(u N (w),c)\d~f N (c) = Q N (z,w) + K , 

“ JC d N 

where K is a universal constant. 

Proof. We showed in EZH p. 779] by an elementary computation that 


IC d N 


log|(ujv(z),c)| log|(MAr(tc),c)|d 77 v(c) = G( \ (u N (z) , U N (w)) \) , 


where 


G(t) = — / e (l Cl l 2 +l c d-)io g | Cl | i 0 g cif + c 2 Vl — t 2 


K 2 Jc 2 


dcidc 2 ■ 


(37) 


(38) 


m 


The computation of the integral (IHH1) was begun in [ SZIl §4], Let us finish it. By (47)-(50) 
SZlj (with A = fr 2 ), we have 


G(e A ) = h + k 2 A + - / log(l - e 2s ) ds 

^ .lr\ 


Since G(0) is hnite, k 2 = 0 and hence 


G(e A ) = h - - / log(l - e 2s ) ds , 


or equivalently, 


Hence, 


G(t) = k 0 - - 


log(l - s) 


ds (0 < t < 1) . 


G(t) = ~[G(t) - k 0 ] . 


7T- 


The lemma follows from EH and (ED with 


„ , . |(^ r iv(-2 ; ), V I / a('11'))| ,, , N -7—r. i 

t = P N (z,w) = r — . rr = (u N (z),U N (w))\ . 


(39) 

(40) 

(41) 

(42) 

□ 
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Proposition 13.21 is an immediate consequence of Lemmas 13.3113.41 □ 

To state the bipotential formula for the volume variance of Theorem II .21 (and number 
variance of Theorem ECU), we let 

* := rT— or - 1 , (43) 

(m — 1)! 

so that Vol 2m — 2 \Zs N n U] = ( Z SN , Xu $)• (If rn = 1 , we set $ = 1.) 

Proposition 3.5. Let ( L,h ) —■> (. M,u ), U c M be as in Theorem \1.?A Then 

Var(Vol 2m _ 2 [Z SN H U]) = — [ B z B w Q N (z,w ) A $(z) A $(w) , 

JdUxdU 

where Qn is given by (ED- In particular, for the case where dim M = 1, we have 


Var (#[Z SN n U})=- [ B z B w Q N (z, w) . 

JdUxdU 


Proof. As in the proof of Lemma [3.31 we let 

y n '■= f-<9<91og|(c,Tjv)|, xu$) = [ -dd\og\(c,ty N )\ A $ = - [ Blog\(c, 'L iV )| A 4> , 

V 77 / Ju 77 71 Jau 

(44) 

where Tjy is given by (I2H1) . As before, we have (Z SN ,xu ®) — y n + NC, where C = 
±f u cvA$ = fVol(U). Hence, 

Var(Vol 2 m _ 2 [Z Sff n U ]) = Var (Z SN , X u $) = Var(V^) . (45) 

By (l44l) (see also the proof of Proposition 12.11) . we then have 


Also by (1441) . we have 


E (Y%)=-[ aioglT^I A<f>. 
71 JdU 


(46) 


E (Tiv) 2 ) = -j [ [ [ [9log |<c, * jv(s)) | A#(z)] [91og|(c,* w (ro))| A*(»)] d-y N (c) . 

77 JC d N JdU JdU 

(47) 

Again following the proof of Lemma 13.31 we write ®n= |'I / v|ma with \un\ = 1, and use (|SS1) 
to decompose (ED into four terms. The hrst term contributes 

—j[ I plog \®n{z)\ a $(z)] [91og|^iv(iu)| A <L(w)] = (EV ^) 2 . (48) 

77 JdU JdU 

The second term vanishes since / log |(c, un{w))\ d^N^c) is independent of w and hence B 
kills it. Similarly, the third term vanishes, since it contains / log |(c, Ujy(z))\ d 7 jv(c), which 
is independent of z. Therefore, 

Var(V^) = [ I I [<91og|(u iV (-),c)| A $(z)] [91og|(u7v(w),c)| A $(w)] ^(c) 

77 Jc d N JdU JdU 

= I B z B w (l log |(ujv(^),c)| log \(u N (w),c)\ $(z) A ®(w) ) <Ixn{c) . 
77 Jc d N \JdUxdU J 

The formula of the proposition then follows from Lemma [3.41 □ 
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4. Off-diagonal asymptotics and estimates for the Szego kernel 


In this section, we use the off-diagonal asymptotics for IIjv(z, w) from ) SZ 2 i to provide the 
off-diagonal estimates for the normalized Szego kernel Pn{z, w) that we need for our variance 
formulas. Our estimates are of two types: ( 1 ) ‘near-diagonal’ asymptotics (Propositions 

Oil for Pn(z,w) where the distance d(z,w ) between z and w satisfies an upper bound 
1 /2 

d(z,w) < b (b G M + ); (2) ‘far-off-diagonal’ asymptotics (Proposition 14.211 where 

d(z,w)>b('-^) 1/2 . 

As discussed in gQ(cf. [zil EZU EZ 2 ] ), we obtain the asymptotics by identifying the line 
bundle Szego kernel IItv with a scalar Szego kernel Iljv(ai, y) on the unit circle bundle X C 
L ” 1 —> M associated to the Hermitian metric h. Given Zo G M, we choose a neighborhood 
U of zq, a local normal coordinate chart p : U,z 0 —> C m , 0 centered at z 0 , and a preferred 
local frame at Zq, which we defined in jSZ2| to be a local frame such that 

l|e i W||» = l-l|W*)f+ ••• . (49) 

For « = («!,..., u m ) G p(U), 6 G (— 7 r, 7 r), we let 

e ie 

p(ui, ■ ■ ■, Um, 0) = | . _ 1( TTT e* L {p-\u)) G X , (50) 

I e L\P \ U )>\h 

so that (mi, ..., u m , 0) G C m x M give local coordinates on X. As in [SZ2j, we write 

n z °(u,e-,v,p) = r i N (p(u,e),p(v,(p)). 

Note that 11^ depends on the choice of coordinates and frame; we shall assume that we are 
given normal coordinates and local frames for each point Zo G M and that these normal 
coordinates and local frames are smooth functions of Zq. 

The scaling asymptotics of 11^ (w, 9] v, <p) lead to the model Heisenberg Szego kernel (f?4l) 
discussed in H 2 . 2 l for the Bargmann-Fock space of functions on C m . We shall use the following 
(near and far) off-diagonal asymptotics from jSZ2| : 


Theorem 4.1. Let ( L,h ) —> (M,u) be as in Theorem [7~^ and let zq G M. 
above notation, 


i) 


tvt—ttitjzo I u 0 . v 

NlpV’ V’ ,/n ’ N> 


Then using the 


= nf {u, 9 ; v, ip) 1 + Eti N-P 2 Vr {u , v) + N~( k+i y 2 R Nk {u, v ) 


where the p r are polynomials in (u,v) of degree <5r (of the same parity as r), and 

\X J R Nk (u,v)\ < C jkeh N e for |u| + \v\ < b^logN , 

for e, b G M + , j, k > 0. Furthermore, the constant Cj k£ b can be chosen independently 
of z 0 . 

ii) For b > \Jj + 2k + 2m, j,k > 0, we have 


v{n N (z,w)\ = o(N~ k ) 


uniformly for d(z,w ) > b 



N 
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Here Vf = (VA )- 7 is the j-th iterated horizontal covariant derivative; see (II 7jl . Theorem 
IP is equivalent to equations (95)-(96) in jSZ2j , where the result was shown to hold for 
almost-complex symplectic manifolds. (The remainder in (i) was given for v = 0, but the 
proof holds without any change for v 7 ^ 0 . Also the statement of the result was divided into 
the two cases where the scaled distance is less or more, respectively, than N 1//6 instead of 
vTog N in the above formulation, which is more useful for our purposes.) A description of 
the polynomials p r in part (i) is given in [S72| . but we only need the k = 0 case in this paper. 
For the benefit of the reader, we give a complete proof of Theorem 14.11 in EH below. 

Remark: The Szego kernel actually satisfies the sharper ‘Agmon decay estimate’ away from 
the diagonal: 

V j U N (z, 9■ w,p) = 0 (e-AjVNd(z,w)^ 5 ■ > o ( 51 ) 

In particular, 

|njv(z,«;)| = O (e- AVWd ( z ’ w ^ . (52) 

A short proof of m is given in jBel Th. 2.5]; similar estimates were established by M. Christ 
EH. H. Delin jDej, and N. LindholmjLi]. (See also |DLML iMMj for off-diagonal exponential 
estimates in a more general setting.) We do not need Agmon estimates for this paper; instead 
Theorem rm suffices. 

We now state our far-off-diagonal decay estimate for P N (z,w), which follows immediately 
from Theorem IP ii) and the fact that n^;?, z) = N m (l + 0(N x )) (by fZij o r Theorem 

EH i))- 

Proposition 4.2. Let ( L,h ) —> (M, u) be as in Theorem II.,4 and let Pn{z,w) be the 
normalized Szego kernel for H°(M, L N ) given by 42(i|) . For b > \Jj + 2k, j, k > 0, we have 


V J Pn{z,w) = 0(N k ) uniformly for d(z,w ) > b 



The normalized Szego kernel P N also satishes Gaussian decay estimates valid very close 
to the diagonal. To give the estimate, we write by abuse of notation, 

P N (z 0 + u,z„ + v) :=P„U>- (u),p- M)- n ^ 0 . Ut0)1/ ^ (Vt0;Vt0)1/ , - 

As an immediate consequence of Theorem EHi), we have: 


Proposition 4.3. Let P n (z,vj) be as in Proposition \f.^\ and let z 0 e M. Forb,e > 0, j > 
0, there is a constant Cj = Cj(e, b), independent of the point z 0 , such that 


M z ° + A' 2 » + 7 f) = + «„(«, t-)] 

\V j R N (u,v)\ < CjN~ 1/2+£ for \u\ + |u| < by/logN . 


As a corollary we have: 


Proposition 4.4. The remainder in Proposition^^ satisfies 

\R n (u,v)\ < ^ \u—v\ 2 N~ 1 / 2+£ 1 |Vi2jv(w)| < C 2 \u—v\ A^ _1//2+£ , for | m | + | t | < by/logN. 
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Proof. Since -P/v (^o + u, zq + v) < 1 = Pn (zq + u, zq + u), we conclude that Rn(u,u) = 
0, d-RAr|( U)U ) = 0, and thus by Proposition 14.hi 

\WRn(u, u)| < sup | V 2 /?at(u, (1 — t)u + tv)\ \u — v\ < C 2 \u — v\N~ l / 2+£ . 
o<t<i 


Similarly, 


Rn(u,v) | < 


- sup | X 2 R N (u, (1 — t)u + tv) | |u — v \ 2 < \u — v \ 2 N l ^ 2+e . 
2 o<t<i 2 


□ 


4.1. Proof of Theorem 14.11 In this section, we sketch the proof of Theorem 14.11 The 
argument is essentially contained in [SZ2j . but we add some details relevant to the estimates 
in Theorem im 

The Szegb kernels IIjv(x, y) are the Fourier coefficients of the total Szego projector II(x, y ) : 
£ 2 (X) —► 7i 2 (X ); i.e. U N (x,y) — f e~ lN 9 H(e l 9 x, y) dd. The estimates for II v(z, w) 
are then based on the Boutet de Monvel-Sjostrand construction of an oscillatory integral 
parametrix for the Szego kernel: 


B (x,y) = S(x,y)+E(x,y) , 

with S(x, y) = J 0 °° e lt ^ x,yS) s{x, y, t)dt , E(x,y) 6 C°°(X x X). 


(53) 


The amplitude has the form s ~ YlkLo^™ ks k{x,y) G S m (X x X x M + ). The phase function 
if is of positive type, and as described in DSZ2j . is given by: 


if(z, 6 ,w,tp) 


a{z,w) ei(0 _ v) 

V a ( z )V a ( w ) 


(54) 


where a G C°°(M x M) is an almost holomorphic extension of the function a(z,z) := a(z) 
on the anti-diagonal A = {(z, z) : z £ M}, i.e., da vanishes to infinite order along A. We 
recall from m that a(z) describes the Hermitian metric on L in our preferred holomorphic 
frame at zo, so by m, we have a(u) = 1 + |u| J + 0 {\u\ ; ), and hence 

a(u, v) = 1 + u ■ v + 0(\u \ 3 + |u| 3 ) . (55) 


For further background and notation on complex Fourier integral operators we refer to [HSZ 2 ] 
and to the original paper of Boutet de Monvel and Sjostrand ra- 

As above, denote the N -th Fourier coefficient of these operators relative to the S 1 action 
by IIjv = Sn + En. Since E is smooth, we have En(x, y) = 0(N~°°), where 0(N~°° ) denotes 
a quantity which is uniformly 0(N~ k ) on X x X for all positive k. Then, En{z,w) trivially 
satisfies the remainder estimates in Theorem o 

Hence it is only necessary to verify that the oscillatory integral 

/»27T S*OQ /»27T 

S N (x,y)= e- iN 9 S(e i 9 x,y)dd= / e- iN 9 +it ^ eiex ’ y) s(e i 9 x,y,t)dddt (56) 
Jo Jo Jo 
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satisfies Theorem rm This follows from an analysis of the stationary phase method and 
remainder estimate for the rescaled parametrix 


U , 0; — 7 =, 0 1 = N 


qzo 

N Wn’^Vn 


f* 00 /*2n 


1 0 JO 




Vb V 

6 ; —=, 0, Nt ) dddt , 


Wn Vn' 


(57) 

where we changed variables t > Nt. For background on the stationary phase method when 
the phase is complex we refer to 11 lol . We are particularly interested in the dependence of 
the stationary phase expansion and remainder estimate on the parameters (u, v ) satisfying 
the constraints in (i)-(ii) of Theorem im 

To clarify the constraints, we recall from jSZ 2 j (95) that the Szego kernel satisfies the 
following far from diagonal estimates: 

AT 1 / 6 


| V^II A r(^, w)| = 0(N K ) for all j, K when d(z,w ) > 


y/N 


Hence we may assume from now on that z = z 0 + -d=, w = z 0 + -^= with 


it 


+ Id < 5N 1/6 


(58) 


(59) 


for a sufficiently small constant S > 0 (see| 

By (Evil) (tool) , the rescaled phase in m has the form: 


T := t'lp 




Wn' W 

and the A r -expansion 


, 0 ) — 6 = it 


1 - 


U V 

Vn 1 Vn 


jo 


U U 

Vn 1 Vn 


V V \ 

Vn' vdv ) 


-e 


(60) 


it 


JO t 4--ni’t U V 


T = it[ 1 - e ] - 9 - —i> 2 (u, v)e l + tR 3 ( 


jo 


'Vn' Vn 


(61) 


where 


i^ 2 (u, v) — u ■ v - (\u\ 2 + 1 17 1 2 ) = —| u — v\ 2 + i Im (u ■ v ) 


is the phase function of After multiplying by iN, we move the last two terms of (ED 
into the amplitude. Indeed, we absorb all of exp-f^ + iNR^)te ld } into the amplitude so 
that ED is an oscillatory integral 

/»oo 

N / 


'0 Jo 


with phase 


and with amplitude 


e iN*w) A fa Q. zQj v y Mt + o(N-°°) 
V(t,9) :=it{l-e iG ) -6 


(62) 

(63) 


A(t,9;z 0 ,u,v) := (64) 

The phase T is independent of the parameters (u,v), satisfies Re (iT) = —t( 1 — cos 6 ) < 0 
and has a unique critical point at {£ = 1 , 6 = 0 } where it vanishes. 
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The factor e teW is of exponential growth in some regions, ffowever, since it is a 
rescaling of a complex phase of positive type, the complex phase iiVT + te ld, ^ 2 (u,v) is of 
positive type, 

Re (iiVT + te* 6 , ' 02 (w, v)) < 0 (65) 

once the cubic remainder Nte ld P^(-^=, -)=) is smaller than iNty + te ld 'ip 2 {u, v ), which occurs 
for all (t , 9, u, v ) when (u, v ) satisfy (I^TJl) with 5 sufficiently small. 

To estimate the joint rate of decay in ( N,u,v ), we follow the stationary phase expansion 
and remainder estimate in Theorem 7.7.5 of [Hoj . with extra attention to the unbounded 
parameter u. 

The first step is to use a smooth partition of unity {pi(t, 6), p 2 (t, 6)} to decompose the 
integral m into a region (1 — e, 1 + e) t x (—e, e)g containing the critical point and one over 
the complementary set containing no critical point. We claim that the p 2 integral is of order 
N~°° and can be neglected. This follows by repeated partial integration as in the standard 
proof together with the fact that the exponential factors in m decay, so that the estimates 
are integrable and uniform in u. 

We then apply jHoj Theorem 7.7.5 to the p\ integral. The first term of the stationary 
phase expansion equals N m e te and the remainder satisfies 

\Rj(P 0 ,u,v,N)\<CN~ m+J su P | D« ePl A(t,9-,P 0 ,u,v)\. ( 66 ) 

H<2J+2 t,d 

From the formula in (EH) and the fact that s is a symbol, A has a polyhomogeneous expansion 
of the form 


A(t,9;P 0 ,u,v) = p^ey^^N* 


K 


N~ n / 2 f n (u, v-1, 9 , P 0 ) + Rk{% v , t, 9) 


n =0 


|V^fl m (u,B)| < C, k€b e ’(W +')»-+). 


(67) 


The exponential remainder factor e £ ^ u l 2+ ^ 2 ) comes from the fact R ee 19 ^ = cos^Re^ — 
sin 6dm 'tp with Re^ < 0 and | sin 9 < e on the support of p±. ffence, the supremum of the 
amplitude in a neighborhood of the stationary phase set (in the support of pi) is bounded 
by e £ l Im b 2 |_ The remainder term is smaller than the main term asymptotically as N —> oo 
as long as (u, v ) satisfies (jHHl) . Part(i) of Theorem 14.11 is an immediate consequence of (ElD 
since e £ ^ 2+ ^ 2 ' > < N £ for |m| + |u| < ydoglV. 

To prove part (ii), we may assume from m-m < |w| + \v\ < 5N 1 / 6 . In 

this range the asymptotics (E3) are valid. We first rewrite the horizontal ^-derivatives 
as Uj derivatives, which for L N have the form \/~N — NAj(^=) and thus V/ t contributes 

a factor of y/N. We thus obtain an asymptotic expansion and remainder for Vj)IIjv(£, w) by 
applying V J h to the expansion (i) with k = 0: 


n?(u, 9] v, if) [l + N 1/2 R N0 (u, u)] . 








NUMBER VARIANCE OF RANDOM ZEROS 


21 


The operator V 3 h contributes a factor of N ^ 2 to each term, and thus 
|V{njv(z,iy)| = O 

= 0(N~ k ) uniformly for \u\ 2 + \v\ 2 > (j + 2k + 2m + e')logN , 

where s' = (j + 2 k + 2 m + l)e. □ 


5. The sharp variance estimate: Proof of Theorem 11,31 

We first give the proof of Theorem II .31 which uses the same method as the proof of 
Theorem o but has simpler computations. 

We begin with some off-diagonal asymptotics for the function Qn = G o P N defined in 
(EH) • By Proposition 14.21 we see that 


Q n (z,vj) < 


1 log(l a) 


4tj-2 


ds = O 


N 2r ' 


5 V log N U;o ^ 

tor d{z,w) > —~^=—, ( 68 ) 


with b > \j 2 m. 

Next we show the near-diagonal estimate 

Qn(z 0 ,z 0 + — 0 (e _ sM 2 ) + 0 (N~ 1/2+£ ) , for |u| < by/logN. (69) 

To verify (jHHll . we apply Proposition 14.41 Since Pn{zo,Zo) = 1 and G'(t) —> oo as t —> 1, we 
need a short argument: let 

Aiv = -log P N . (70) 

Recalling (1251) , we write, 

2 P°° 


F(X) : = G(e~ x ) = 


2 tt 2 


log(l — e ) ds (A > 0) , 


(71) 


so that 


Qn — F o \ N . 


By Proposition 14.41 


Aw 


^0, Z 0 + 



1 

2 


u|- + R n (v) , 


where 

Rn — — log(l + Rn) = OdvfN- 1 ' 2 **) 
By (ED and Proposition 14.41 


(72) 


l^(A)| = 


2 vr 2 


log(l 


-2Aa 


< 


2vr 2 


max 


log A’ 1 


= 01 + log + — 
\v 


(73) 
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Since i|w| 2 + Rn(v) = \v\ 2 (| + o(iV)), it follows from (1721) - (I7H that 


Qn(zq,Zq + ~^=j — F ^l^l 2 Rn( v ) 


2' 

= F[±\v\ 2 )+0 


1 + log H 


Rn(v) 


= G(e"3H 2 ) + 0(iV- 1 / 2+£ ) , for \v\ < b^logN, 

which gives flfihjl . (This computation also shows that Qn is C 1 and has vanishing first 
derivatives on the diagonal in M x M.) 

Next, we note that G(t) = -A (t 2 + £ + £ + ••• + , and hence 


- it l 


k =1 


E 


D -k\v\ 2 

k 2 


■ dv 


7r 


7r 


m— 2 


47T 2 ^ k m+2 

fc=l 


C(m + 2) 


By Proposition El we have 


Var(E Sjv , </?) = / l(z)iddip{z) , 


i M 


where 


We let 


liz) = 


' {z}xM 


Qn{z, u>) iddip(w) . 


(74) 

(75) 

(76) 


n = —u r 
ml 


denote the volume form of M, and we write 

iddip — 7 /; Q , if; e C 1 . 


(77) 


To evaluate Z(z 0 ) at a hxed point zq E M, we choose a normal coordinate chart centered at 
Zq as in (gj By (jHHI) and (17(1 - (1771) . 


Az 0 ) = 


l\v\<b^/logN 


Qn ( Zq, Zq + 


y/N 


fp\Zo + 


y/N 


( Zq + 


y/N 


+ 0 


N 2r ' 


(78) 


We recall that u> = log a = | dd [|z | 2 + 0(|z|'' : )] in normal coordinates. Hence 


H ( Zq + 


An) ml [2N 

for |v| < by/logN, where 


^-dd\v \ 2 + 0 


_HA 

AT3/2 J 


N r 


£Ie(v) + O 


log n\ 
— 


(79) 


= T ('zdd\v\ 
ml \ 2 


nn rn 


| J -dvj A dvj 


3 = 1 
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denotes the Euclidean volume form. Since <p E C 3 and hence ip(z + ^f=) = i/j(z) + 0 (\v\/VN), 
we then have by © and CED-m 
1 


A z o) = 


N' n \_J |u|<fe v / Iog _ /V 

x{^(n)+0(iV- 1 / 2+£ )} 


\G(e~A\ 2 ) + QiN- 1 ^)} {A z o) + 0(N~ l ' 2+£ )} 


+ 0 


N 2r 


Since G(e ^M 2 ) e £} by (I7i|) . we have 

V’(^o) 


A Z o) = 


N r ‘ 


-J\v\<by/log N 


G(e“5 M )n E (v) + O (iV 1 / 2 +£ (logiV) m ) 


Since G(e A ) = 0(e 2A ) and hence 


G(e~Al 2 )n E ( v ) = 0(N~ 2m ) , 


'\v\>by/log N 


(80) 


(81) 


(82) 


we can replace the integral over the ( 6 -^/log7V)-ball with one over all of C m , and therefore 


A*o) = 


A z q) 

N m 


G(e~A\ 2 ) n E (v) + 0(N~ l / 2+£ ') 


_Jc n 


A z o) 

N m 


Hi r , 


+ 0(N~ 1 / 2+£ ')\ , (83) 


for all s' > £, where K m = T ' r '\ Q{rn + 2) by (jTjj) . Therefore, by m and 


Var (Z 8 ff ,(p) = 


N Tl 


' M L 


Km + 0(N~ 1 l M )} 1>(Z) 2 Q(Z) . 


(84) 


Since 


dHU) yields the 


A z ) 2 A z ) = / \dd<p\ 2 Q = \\dd<p\\l, 


' M 


variance 


formula of Theorem o 


□ 


6. Variance of zeros in a domain: Proof of Theorems 11.1 1 PH 

Following the approach of we now prove Theorem II .21 and, as a consequence, we also 
obtain Theorem o which is the one-dimensional case of Theorem PI 
By Proposition 13.51 we have 

Var(Vol 2 rn . 2 [Z SN n U]) = f T A <f> , (85) 

JdU 

where <f> = — Arzuj m ~ l and 

(m—1)! 

T(z) = -B z [ B w Q n (z,w ) A A w ) ■ ( 86 ) 

J{z}xdU 

From (J2H) and Proposition 14.21 we conclude that 

B z B w Q N (z, w) = 0(N~ m ) , for d(z, w) > b\f^- , (87) 

where we choose b = \j2m + 3. Thus, we only need to integrate over a small ball about 
z in dU when using to compute the variance. 
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To evaluate T(zo) at a fixed point z o G dU, we choose normal holomorphic coordinates 
{zj} at zo, defined in a neighborhood V Of zq. By dHSD-dHZD, we have 


T A 0 = -Y.il 

i* V s 


d 2 


au dzjdwk 


Qn{z,w) 


dwk A <&(w) ] dzj 


Z = ZQ 


E 

j,k 


d 2 


O 


[ Z0+w edU:\w\<b^§^} dzjdwk 

1 

iW 


— Qn(z 0 + z,z 0 + w) 


dwk A $(w) ] dzj 


2=0 


( 88 ) 


As in the proof of Proposition 11.81 we write Qn = F o An, where A n = — log Pat and F 
is given by (I7H) . By Propositions 14.1114.41 


N ~ 

A n (z 0 + z,z 0 + w) = —\z - w\ 2 + R n {z, w) , 


where Rn(z, w) — — log 1 + Rn (y/N z, '/N 


w 


satisfies the estimates 


By (ED, 


| J Rj V (z,iu)| = 0{\z-w\ 2 N l / 2+e ), \VR n (z,w)\ = 0(\z — tc|A rl / 2+e ), 
|V 2 Pjv(-,w)| = 0(N l / 2+e ), for \z\ + M < 


f '( A) = Pio g( i-e-), m = 


We now let A = An(zo + z, Zo + w). By (IHIl - flHUl) . 

\F'(X)\ < 2^2 max flog1^ = O (l + log 


A' 


z — w 


F"( A) < 


= O 


2ir 2 \ \iV|^ — w\ 

Hence, for \w\ < b\J— jfX, we have 


for \z\ + \w\ < b 


log AT 


N 


d 2 

0 _ 0 _ Qn{Zq + z, z 0 + w) 

OZjOWk 


z =0 




2 = 0 


--Nwj + 0(\w\NF 2 + £ ) 


= F"( A) 

+F'(A) ■ 0(N 1/2 +‘ 


-Nw k + 0{\w\NF 2 +z 


- ^N 2 F , \X)w j w k + 0(NF 2 + £ )(l + log + -^] . 


Furthermore, since 


-F< 3 )(() = -csch 2 (<r 2 , 


(89) 


(90) 


(91) 


(92) 
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we have 


and hence 


F"{ A) = 


n 2 (e N \ z w \ 2 — 1 ) 


+ 0{\z-w\~ 2 N~ 3/2 + s ) , 


d 2 


a- Qn{z o + Z, M) + w) 

UZjUWk 


z =0 


-AT 2 


rtUjlUfc + 0(iV 1/2+£ ) ( 1 + log + — 
' ' \w\ 


(93) 


4 7r 2( e A|«)| 2 _ 

We note that under our hypothesis that dU is piecewise C 2 without cusps, we have the 
estimate 

1 
2 


J {zo+w£dU:\w\<5} 

Substituting into dHHD, we obtain 


(l + log + dVol(w) = O (5 2m_1 I log<51) , for 6 < 


(94) 


™- 5 (S L 


w.jWkdwk 


{zo+wedU:\w\<by/ 1 -^} e N \ W \ 2 - 1 


/\$(w) + 0(N 


—m+l+2e 


) dzj. (95) 


We first consider the case where dU is C 2 smooth (without corners). We can choose our 
normal coordinates {zj} about z 0 so that the real hyperplane (Im z { = 0} is tangent to dU 
at Zo = 0. We can then write (after shrinking V if necessary), 

U fl V — {z e V : Im z\ + <p(z) > 0} , 

where tp : V —> M is a C 2 function of (Rezi, z 2 ,..., z m ) such that 99 ( 0 ) = 0, d(p( 0) = 0. 

We let 

r(v) = (vi + icp(v),v 2 ,..., Vm) , 

so that dU = (Im V\ = 0} in terms of the (non-holomorphic) v coordinates. We make the 
change of variables 


in the integral (USD: 


N 2 
47T 2 


tw = £ 

j,k 

where 

i6lx C m_1 : \v\ <{b- 


Bi 


w = tn(v) := t 

WjWkdv)k 
e N M 2 - 1 


T, 


N 


Vn 


A <f)(tc) 


+ O (W m+1+2e ) dzj, 


(96) 


c B 1 c <! v e M x 


C 71 - 1 : M < (6 + l)v/logiv). 


To evaluate the integrand in (Infill , we first note that 


w 1 = 


t’i 

v/fV 


iV 


+ O —— , dw 1 = —= dv\ +0 — , M7 2 = 




Thus iV|-u ;| 2 = |n| 2 + O ( ), and hence by 


e N\w\ 2 _ e \v\ 2 _ 


N 


O [ ^ ) 0 (\v\- 2 ) = 


V2 

y/N' 


• • 1 


y/N 


e H 2 - 1 


+ 0(N ~ 1+£ ), 
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for |u| < 2 b^/logN. Finally, we have 


AT3/2 J 


m— 1 


N- 


'771—1 


^W + oh/^j 


on S 2 ™" 1 , where 


®e(v) = 


(m- 1)! \2 


i 


<9<9|u| 


771—1 


Therefore, (1HHI) becomes 

m 

T{z 0 ) = jV-™+ 3 / 2 ^ 

i=i 

We note that 


'Bi 


-i e.H 2 - 1 


1 / ^ 

47T 2 J d 2m-l eH 2 — 1 


<9|u| 2 A ^(u) = 


<9|n| 2 A<h B (n)+ 0(iV 1 / 2+2£ ) 


-i e.H 2 - 1 


^Vol MxC m-i(n) . 


Since 


m 


, , |2 do: = 0(AT- b2+1 ) , 

^xgR n :|x|>f)yiogN' ^ ^ 

we can replace the integral in m with an affine integral, and hence 


(97) 


dzj . 


(98) 


(99) 


m 1 /■ 

T(«b) = ^“ m+3/2 E — / 

pT L 47r2 A 


where 


V 3 V 1 

RxC m_1 — 1 


dVol RxC _ 1 (n) + 0(iV- 1 / 2+2£ ) 


= AT- m +3/2 ^ ^ + 0 ^-m+ 1+ 2 £ ) ^ 


XA 




4(/r 2 


»-! elb 2 - 1 


dx = 


\x\ 


2-^171-1/2 roo r 2m 


47t 2 (2 m — 1) r(m — 1/2) J 0 e r — 1 

m—5/2 °° /■oo 

4F(m+ 1/2) 

5/2 


47r 2 (2m — 1) j K 2m-i e^l 2 — 1 
dr 


dx 


^ /»oo 

£ / e~ kr \ 2m dr 

k =l ' /o 

f- r(m+l/ 2 ) = ^ m - 5/2 7 1\ 

4T(m+ 1/2) “ 2fc™ +1 /2 8 H 2/ 


dz,- 


( 100 ) 


Substituting w into dHSD, we obtain the formula of Theorem PI which completes the 
proof for the case where dU is smooth. 

We now consider the general case where dU is piecewise smooth (without cusps). Let S 
denote the set of singular points (‘corners’) of dU, and let S/v be the small neighborhood of 
S given by 

b'y/]ogN\ 

—JlT J ’ 

where b' > 0 is to be chosen below. We shall show that: 


Sn — s z E dU : dist(z, S) < 
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i) m holds uniformly for zq G dU \ Sn] 

11) sup \T(z)\=0(N- m+3 / 2+£ ) . 

z&dU^S 

Since Vol 2m _iSjv = O ( , the estimate (ii) implies that the integral in 425D over the 

small set is negligible and hence 

Var(Vol 2m . 2 [Z SN n U}) = f TAH 0(N~ m+1+2£ ) . (101) 

JdU^S N 

It then follows from (i) and cmD that 


Var(Vol 2m _ 2 [Z Sjv n U)) 


jy—m+ 3/2 
jy-m+3/2 


I'm Vol 
I'm Vol 


2m— 1 

2m-1 


{dU ^ S N ) + 0(N~^ +2e ) 
(dU) + 0{N~^ +2e )] , 


which is our desired formula. 

It remains to prove (i)—(ii). To verify (ii), for each point Zq G dU \ S', we choose holomor- 
phic coordinates {zj} and non-holomorphic coordinates {vj} as above. We can choose these 
coordinates on a geodesic ball V ZQ about zq of a fixed radius R > 0 independent of the point 
Zq, but if z 0 is near a corner, dU will coincide with {Irn V\ = 0} only in a small neighborhood 
of zq. To be precise, we let K Zo denote the connected component of I4 0 DdU \ S containing 
Zq. Then we choose ip G C 2 ( I4 0 ) with <p( 0) = 0, dip( 0) = 0, such that 


{z G V : Im Z\ + <p(z) = 0} D K Zo . 


( 102 ) 


Choose N 0 > 0 such that b' < R ; then 

iwedU: d(z 0y w) < j C W 0 , for N > N 0 . 

Then for N > Nq, the integrals © and flSHI) hold, except now they are over a piecewise 
smooth hypersurface B N of the (6'v / logiV)-ball in C m instead of the linear hypersurface 
B 2 ™-\ Since 

Vjd\v \ 2 ^ \v \ 2 ^ 

eH 2 - 1 “ eH 2 - 1 “ ’ 

where the above norm is respect to the Euclidean metric on T*( (C^), it follows from (ITTH1) 
(with Bn replaced by Bn) that 

ITM < m +^° (1) jy-m+3/2 Vo^_ 1 (B Jf ) , (103) 

where Vol' 6 ' denotes Euclidean volume. Since dU is piecewise smooth, we see that 
Vol^S*) < N™- 1 ' 2 [l + o(l)} Vol 2m _ 1 {wedU:d(z 0 ,w)<b'^^} 

= N™” 1 ' 2 O ([^] m “ 1/2 ) = O ((logiV)” 1-1 / 2 ) . 

Combining m-m, we obtain the bound (ii). 


(104) 
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To verify (i), we let 


C = sup 


dist(z, S) 


zedu^s dist (z,dU \ K z ) ' 

We recall that our assumption that dU ‘has no cusps’ means that U is locally C 2 diffeomorphic 
to a polyhedral cone, which implies that C < +oo. We now let b' = Cb, where b = \/2m + 3 
as before. 

Consider any point z o G dU \ Sn, N > N 0 . Then 


dist(z„, BU x KJ > dist ^' S) > . 

1 c - cVN Vn 

Thus by our far-off-diagonal decay estimate (ED, the points in dU \ K Zo contribute negligibly 
to the integral in (19 til) , so that integral can be taken over ( K zo ), or over the linear 

(b^log N)-ball B n . Then (HOOD follows as before. 

Thus we have verified (i)-(ii), which completes the proof of Theorem 11.21 for the general 
case where dU has corners. □ 


7. Asymptotic normality: Proof of Theorem 11.41 

The proof is a combination of Propositions 14.21 and 14.41 with a general result of Sodin- 
Tsirelson jftTj on asymptotic normality of nonlinear functionals of Gaussian processes. Fol¬ 
lowing eh. we define a normalized complex Gaussian process to be a random function w(t) 
on a measure space (T, /i) of the form 

w{t) = ^(t) , 

where the Cj are i.i.d. complex Gaussian random variables (of mean 0, variance 1), and the 
fj are (fixed) complex-valued measurable functions such that 

X] \fM 2 = 1 for a]1 t e T - 

We let wi,w 2 ,w 3 , ... be a sequence of normalized complex Gaussian processes on a finite 
measure space (T, p). Let f(r) G £ 2 (K + ,e -r ! 2 rdr ) and let -0 : T —> M be bounded measur¬ 
able. We write 

Zfj(w N ) = [ f(\w N (t)\)i/;{t)dfji{t). 


Theorem 7.1. jSTl Theorem 2.2] Let p^(s,t) be the covariance functions for the Gaussian 
processes w^ft) ■ Suppose that 

. f T f T \p N (s,t)\ 2a 'i/j(s)fj(t)dfi(s)dfi(t) 
jvi 1 ^ snp seT J T \p N ( s ,t)\dn(t) 

for a = 1 if f is monotonically increasing, or for all a G Z + otherwise; 


n 


lim sup / |pjv(s, t)\dfx(t) — 0. 


N—>oo 


seT JT 


Then the distributions of the random variables 


Av 


EZ* 


\/TTzt) 













NUMBER VARIANCE OF RANDOM ZEROS 


29 


converge weakly to J\f( 0,1) as N —> oo. 


We apply this result in the case fir ) 
Gaussian processes 


w N (z) 


= logr and (T, /i) 
s N (z) 


(M, 0), with the normalized 


where Stv is a random holomorphic section in H°(Ad, L N ) with respect to its Hermitian 
Gaussian measure. The covariance kernel of this Gaussian process is Pn(z,w). Further, we 
let (p be a C 3 real (2 m — 2)-form and we write ddip = i/jQ as before (and hence ij) G C 1 ), so 
that 

Z%{w N ) = (Z SN ,<p) = I log \sN(z)\ddlf(z) 

J M 

is the smooth linear statistic of integration of the fixed test form (p over the random zero 
set. (This was the application of interest in jSTj . where they considered random functions 
on C, CP 1 , and the disk.) 

By Proposition 12.11 we have 


E Z% = ^ dd log njv(«, z) + Nu, 

Z7T 

hence 

Z*(w n )~EZ*= f logb== gS^ 

Jm ^yYn\z.,z) 

To apply the theorem it suffices to check that Pn(z,w) satisfies conditions (i)-(ii). We 
start with (ii): by Proposition 14.21 


lim sup / 

N ^°° Z£M J d {z,w)>y/^°Ki 


\P N (z,w)\dV 0J (w ) = 0. 


On the other hand, since \Pn(z,w)\ < 1 it is obvious that the same limit holds for d(z,w ) < 


blog A 

N ■ 


To check (i), we again break up the integral into the near diagonal d(z } w) < and 

the off-diagonal d(z,w ) > \J bl °^ N ■ As before, the integrals over the off-diagonal set tend to 
zero rapidly and can be ignored in both the numerator and denominator. 

On the near diagonal, we replace P/y by its asymptotics in Proposition 14.31 The asymp¬ 
totics are constant in z and with uniform remainders, so the condition becomes 


lim inf 

N—>oo 


full 


|«|< v/blog N 


e |u|2 |l + R n (u)\ 2 ^(z + -^)'ijj(z) dudVu(z 


h 


|n|< v / 5Tog7V 


e 2 M 2 |1 + R n (u)\ du 


> 0. 


Since ifj G C 1 , the ratio clearly tends to 2 rn j M ij){z) 2 dV u > 0, completing the proof. 


□ 
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